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ABSTRACT 

We discuss the nature of the velocity dispersion vs. size relation for molecular 
clouds. In particular, we add to previous observational results showing that the veloc- 
ity dispersions in molecular clouds and cores are not purely functions of spatial scale 
but involve surface gas densities as well. We emphasize that hydrodynamic turbulence 
is required to produce the first condensations in the progenitor medium. However, as 
the cloud is forming, it also becomes bound, and gravitational accelerations dominate 
the motions. Energy conservation in this case implies \E g \ ~ Ek, in agreement with 
observational data, and providing an interpretation for two recent observational re- 
sults: the scatter in the Sv — R plane, and the dependence of the velocity dispersion 
on the surface density Sv 2 /R oc E. We argue that the observational data are consistent 
with molecular clouds in a state of hierarchical gravitational collapse, i.e., developing 
local centers of collapse throughout the whole cloud while the cloud itself is collaps- 
ing, and making equilibrium unnecessary at all stages prior to the formation of actual 
stars. Finally, we discuss how this mechanism need not be in conflict with the observed 
star formation rate. 

Key words: ISM: general - clouds - kinematics and dynamics - turbulence - stars: 
formation 



1 INTRODUCTION 

Almost 30 years ago, ILarsonl (|l98ll) suggested that two scal- 
ing relations exist for molecular clouds, one for the velocity 
dispersion Sv and the other for the mean density p. These 
have the form 

pcxiT, (1) 
Sv oc R p , (2) 



where R is the size of the cloud. The most commonly ac- 
cepted values of the exponents arc] a ~ — 1 and P ~ 0.5 (see 
ISolomon et all Il987l: iBlit j ll993l : lBallesteros-Paredes et all 
120071 : iMcKee fc Ostrike3l2007l . and references therein). The 
first relation implies that the mean column density of the 
gas, E = pR, is roughly constant for the whole ensemble 
of clouds. This re sult has been challenged by various au- 
thors; for instance, iKegell (|l989h point ed out that i t may be 
the result of various selection effects. IScalol (|l990h showed 
that the study by ISolomon et al] (|l987l ) was sensitive only 
to a limited dynamical range of column densities. More- 
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1 Note that the original values for the exponents reported by 
ILarsonl (Il98lh are a ~ -1.1 and f3 ~ 0.38. 
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over, IVazguez-Semadeni et all (|l997T l used numerical sim- 
ulations of a 1-kpc 2 piece of the Galaxy to argue that the 
mean density-size re lation does not hold when no detectabil - 
ity limitations exist. IBallesteros-Paredes fc Mac Low! (|2002l ) 
used three-dimensional simulations to confirm that clouds 
can have different mean densities, and that the density-size 
relation appears when clouds are observed with particular 
tracers, since they must have at least a minimum column 
density in order to be detectable. 

On the other hand, the velocity-size relation (equation 
[2]| has often been assumed to be real. One reason for this 
is probably that the relation is similar to what that which 
might be expected from studies of fluid turbulence: the 
expected value for incompressible turbulence is /3 ~ 1/3, 
while a turbul ent fluid dominated by shocks might exhibit 
J3 ~ 1/2 (e.g., lElmegreen fc Scaloll2004 iMcKee fc Ostrikeil 
120071 . and references therein). As the interstellar medium 
is highly compressible, and subject to strong shocks of stel- 
lar winds, supernovae, spiral arms, etc., one might then ex- 
pect that the velocity dispersion-size scaling relation with 
/S = l/2 should be valid. And of course, supersonic turbu- 
lent simulations, as well as analytical calculations of molec- 
ular clouds performed over the last decade, have supported 
the ide a of a velocity dispersion-size r elation with /3 ~ 1/2 
(e.g., IVazguez-Semadeni et aL I Il997l; iPadoan fc Nordlundl 



1999l;lBallesteros-Paredes. Vazquez-Semadeni fc Scaloffl999l; 
Ballcstcros-Parcdes & Mac Low 2002; Padoan & Nordlund 



2002l ; lKrumholz fc M cKee 2001; iField et alJl200St etc) 



However, i f Larso n's original dataset, as well as that of 
ISolomon etaH (|l987T l were limited in the range of column 
densities observed, why should only one of these correla- 
tions be affected but not the other one? In principle, one 
may expect that the velocity dispersion-size relation is also 
implicitly a re sult for a lim i ted ra nge of column densities. 

Recently, iHever et al.l (|20091 ) observed an ensemble of 
molecular clouds with the 14 m FCRAO telescope. Because 
these data have improved sensitivity and better spectral and 
angular resolutio n, they r each m uch larger surface densitie s 
than the data in ILarsonl i|l98ll ) and ISolomon et all (|l987h . 
IHever et al.l (|20091 ) found that the velocity dispersion does 
not depend simply on size scale, but on the square root of 
the column density as well. They went on to point out that 
the revised relation, 



5v oc Y, L "R 



1/2 pl/2 



(3) 



was that to be expected for clouds in gravitational equilib- 
rium. In addition, while the derived masses for the clouds 
were a factor of a few low er than expected for Virial equilib- 
rium, IHever et all |2009) stated that uncertainties in their 
mass e stimates still a llowed for consistency with equilibrium 
states. ILarsonl (|l98ll ) actually came to a similar conclusion 
even with his limited dataset; he argued that the clouds 
are mostly gravitationally bound and in approximate Virial 
equilibrium. However, since the clouds' column densities 
vary by over three orders of magnitude (from ~ 10 21 cm -2 
in the most diffuse clouds, to ~ 10 24 cm -2 for infrared dark 
clouds), then the simple form of the velocity-size relation, 
equation may not be the most appropriate form for 
molecular clouds. 



To further explore the correlation between the velocity 
dispersion and size for molecular cloud cores, in iJHwe com- 
pile recent data from the literature in order to show that a 
unique trend for the scaling of the velocity dispersion with 
size does not appear to exist, specifically when we include 
recent sensitive observations of massive cores. Instead, rela- 
tion ([3| seems to hold in all cases. 

The question is then, what is the origin of the 
velocity-surface density-radius correlation? In the past, 
it has been proposed that this relationship arises from 
a condition of hydrostati c equilibrium a p plying to the 
clouds and dense clum ps jElmegreenl 1 19891 ; IMcKee fc Tanl 
120031 ; IField et al.l I2010T I . In this case, it is assumed that 
the clouds are confined by an external bounding pres- 
sure, which is estimated from the clouds' column density, 
and that the role of the turbulent motions within the 
clo uds is to provide support. However, as summarized 
by IBallesteros-Paredes. Vazquez-Semadeni fc Scald (|l999T l 
and IBallesteros-Paredes! l|200eJ l. the complex non-linear, 
large-scale, and anisotropic nature of turbulent motions 
implies that they do not necessarily provide support, but 
rather cause continuous morphing and reshaping of molec- 
ular clouds, and contribute, or pe rhaps are even driven by, 
the c louds' gravitational collapse (|Vazouez-Semadeni et all 
2008). In particular, it is difficult to see how such an 
irregular and locally anisotropic velocity field could "know" 
how to adjust the magnitude and orientation of the 
turbulent motions to maintain clouds in approximate 
equilibrium for several free-fall times. Indeed, it is difficult 
to argue that global equilibrium is maintained, given 
evidence for age spreads in stellar populations that are 
smaller than the lateral crossing times of the c louds 
(Ballcstcros-Parcdes. Hartmann & Vazqucz-Scmadcn: 



adena 

Il999l ; lHartmann, Ballesteros-Paredes fc Berginl I2001T ). As 
we discuss in the present contribution, the observed corre- 
lations can be explained as long as the velocity dispersions 
result from gravitational acceleration, without requiring 
equilibrium at any stage prior to the formation of an actual 
star. 



2 VELOCITY DISPERSION VS SIZE 
RELATION? 

Although it is frequently argued that molecular clouds and 
their cor es usually exhibi t a re lation like eq. after the 
study b y Casern fc Myersl 1119951 ). and probably more clearly 
that bv lPlume et al.l (1 19971 ). it became somehow recognized 
that massive cores may exhibit a shallower slope than the 
frequently quoted Sv oc r 1//2 . Moreover, the various available 
datasets have not been plotted all together. 

In Fig. [T] we plot the velocity dispersions as a functio n 
of size for th e dense c ores given bvlCaselli fc Myersl Jl995h; 
iPlume et all (Il997l1; IShirlev et all (|2003T l; iGibson et all 
(2009) and Wu et all 1 20101 1. For comparison, we have in- 
cluded also the data points of the original work of Larson! 
(|l98ll ). as well as the recent data bv lHever et al.l (|2009T ). nei- 
ther of which focused particularly on dense massive cores. 
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The dotted lines in this figure represent Larson's fit to his 
, i.e 



km s 



= 1.11 ^ 

pc 



(4) 



We observe th at, while in g e neral terms the typical CO 
clouds observed bv lHever et al. l|2009l ) lie close to Larson's 
relation, this is clearly not the case for the dense, massive 
cores, which exhibit large velocity dispersions for their rel- 
atively small sizes. Although the deviation is only marg inal 
for the Orion cores observed by iCaselli fc Myers! (|l995l ). it 
is clearer for the more massive compact cores cores reported 
by the rest of the datasets. 

In Fig. [2]we plot the "Heyer relation", Sv/r 1 ^ 2 vs. sur- 
face density E. Unfortunately, the observations of massive 
dense cores, particularly from infrared dark clouds, are very 
recent, and only few cores have independent mass estima- 
tions. From t h e list of references given above, only that by 
iLarsonl (|l98ll ); iHever et all |2009l) and IGibson et al.1 (120091 ) 
list the masses of their observed cores independent of the 
Virial mass, and thus, only these data can b e plotted. In 
this fi gure, the data point s from the samples bv lHever et al.l 
l|2009h and iGibson et ail (|2009l ) are denoted by H and G, 
respectively. The long-dashed and dotted lines respectively 
represent the loci of structures in Virial equilibrium and 
structures undergoing free-f all. From this figure , it is clear 
that the massive cores from IGibson et all (2009) span over 
two orders of magnitude in column density, from ~ 100 to 
a few times 10 4 Mq pc -2 , the latter being a factor of 10 
larger than the maximum column density in the lHever et al.l 
l|2009l ) data. We note that The entire dataset is seen to be 
reasonably well fitted by a relation of the form proposed by 
Heyer and in fact, at face value, seems to agree better with 
the free-fall regime than with Virial equilibrium. 

It is important to note that the free-fall relation actu- 
ally implies larger velocity dispersions than the Virial equi- 
librium one, contrary to the very common interpretation 
that velocities higher than Virial imply that the clouds are 
unbound. From this discussion, we see that they can mean 
chaotic infalling motions instead. 



3 DISCUSSION 

3.1 The Sv — R relationship, a consequence of the 
role of gravity in the formation and evolution 
of molecular clouds 

In order to understand the nature of the 8v — R relation it is 
important to understand how molecular clouds are formed 
and how they evolve. 

In the last few years, a number of au- 
thors have supported the idea that molecular 
clouds are formed out of atomic gas when large- 
scale streams collide at transonic speedfl (e.g., 



2 These ideas have been developed specifically for understand- 
ing star formation in the solar neighborhood, where most of the 
gas is atomic. In other regions such as the molecular ring, large- 



i Ballesteros-Paredes. Vazquez-Semadeni fc Scald 1 19991 ; 
i Ballesteros-Paredes. Hartmann fc Vazquez-Semadenil 
Il999| y The collision nonline ar ly triggers thermal insta- 
bility in the post-shock gas jHennebelle fc Peraultl 1 19991 ) 
with the result that the gas co ols rapidly, producing 
a dense, cold, tur b ulent cl oud l| Walder fc Folinil 199S; 



Kovama fc Inutsukal 2 0021: 



Heitsch et al 



20051 . 



2006; 



Vazqu cz-Scmadcni et al. 2006) which soon becomes Jeans 



unsta ble (|Vazquez-Semadeni et all 120071 ; iHeitsch et al.l 
2008). Since the interstellar medium, rather than homo- 
geneous, is highly structured, such shocke d, cold, dense 
medi um will naturally produce clumps (|Bonnell et al.l 
2006). As the whole cloud collapses, the column density 
increases rapidly, allowing the formation of molecula r 
gas (e.g., IHeitsch fc Hartmannl [20081 ; I Glover et al.l 120101 ). 
From this point of view, molecular clouds are in a global 
state of collapse, with an internal distribution of free-fall 
times due to the fluctuations in the density field i n duced 
by the initial turbulence jBurkert fc Hartman nl 12004 



I Vazquez-Semadeni et al.ll2006l ; Heitsch fc Hartmannl 12008). 
Note that this does not necessarily mean that the entire 
cloud is collapsing along all dimensio ns; for example, in the 
toy m odel of the Orion A cloud by lHartmann fc Burkertl 
l|2007l ). the angular momentum of the cloud prevents 
collapse along the long dimension. 

In this scenario, supersonic motions in molecular clouds 
are driven by gravitational energy as the clouds pro- 
ceed to collapse, a situation that has been supported by 
IHeitsch et alT (2009). who show that the synthetic line pro- 
files observed in numerical simulations of the process resem- 
ble those in observations, e.g., the magnitude of the velocity 
dispersion of 13 CO line profiles, or the (small) core-to-core 
velocity dispersion. 

The key point, h oweve r, was emphasized by 



IVazquez-Semadeni et all (|2007l ). who showed that the 
kinetic and gravitational energies of the collapsing cloud 
develop a Virial-like relationship in complex gravitationally- 
collapsing clouds, in which the absolute value of the kinetic 
energy is very close to that of the gravit ational energy (see 
Fig. 8 in I Vazquez-Semadeni et al.l |2007| )PI In other words, 
the "Virial" relation between kinetic and gravitational 
energy, rather than being an indicator of true Virial equi- 
librium, simply shows the importance of gravity in driving 
muc h of the non-therm al motions. 

IField et al.l (2008) have advanced a model of a "gravi- 
tational cascade" , which essentially captures the mechanism 
observed in the simulations. Their model is analogous to a 
turbulent cascade, but in which the quantity being cascaded 
is mass rather than kinetic energy, and the main driver 
is gravity at all scales. In their model, these authors pro- 
pose that the contracting motions are somehow random- 



scale flows in molecular gas can also produce dense star-forming 
regions. 

3 Note that in Vazqucz-Scmadcn i et al. I fc007l) , the absolute 
value of the gravitational energy is a few times that of the kinetic 
energy. However, in that paper, the latter energy was computed 
exclusively for the dense gas, while, for simplicity, the gravita- 
tional energy was computed for the whole numerical box. 
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Figur e 1. Velocity dispers i on-size space for obs e rved cores. Data co mpil ed from | L arson| ljl98lh ; ICaselli fc Mversl l|l995h ; IPIume et al.l 
(110971 "); IShirlev et alj \200$)- iGibson et al.l (120091 ); iHever et alj d200^ and lWu et al.1 j201ttl . Note that, while large CO clouds from the 
survey by I Hever*^taL exhibit the typical Larson relationship, dense massive cores are located above this relationship. All together 

(lower-right panel), the data do not exhibit a relation between the velocity dispersion and size. 
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Figure 2. Heyer's relation (Sv / r 1 / 2 v s . surface dens i ty S) for the 



clouds reported in lHever et al.1 J2009h ; [Gibson et al.1 1 20091). Note 
that t he massive cores follow the tendency showed by Heyer et al.1 
(2009), but they are at the large column density end of that re- 
lationship, ranging from 100 to ~ 25,000 Mqpc - 2 . 



ized, so that the kinetic energy released by the collapse 
is converted into quasi-isotropic motions, and Virial quasi- 
equilibrium can be established at every scale (see their §3). 
Such Virial equilibrium states require th at the mass frag- 
ments are bound by an external press ure (|Elmegreenlll989l : 
iMcKee fc fanll2003l ; iField et al.ll2010t ). In the latter paper, 
the authors write the Virial equilibrium equation for a non- 
magnetized cloud of gas of mass M and radius R in terms 
of its mass surface density E = M/R 2 as 



on 



GE 



47T 



(•») 



where the first term in the right hand side corresponds to the 
gravitational energy, and the second one is due to the exter- 
nal pressure acting over the surface of t he cloud. However , 
the successive virializations suggested bv lField et al.l (|200St ) 
wi ll not occur if c l oud l ifetimes are too short, as pointed out 
bv lBonnell et al.1 (120061 ). 

More importantly, a cloud undergoing collapse does not 
need to be confined by an external pressure. Instead, the in- 
ternal pressure increases together with the density, and once 
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the collapse has advanced sufficiently, the external pressure 
becomes negligible. Indeed, as mentioned above, these re- 
gions within molecular clouds are known to hav e much large r 
thermal pressures than the ISM mean (e.g., iBlita 11993). 
Thus, rather than Virial equilibrium, the relevant principle 
here is that the total energy of the system (E s + E^) is 
conserved. In this case, the ratio dv 2 /R is given by 



ov 
~R 



2GE. 



(6) 



Moreover, if the collapsing scenario applies to all scales 
within molecular clouds, we expect eq. ([6| to be valid not 
only for massive cores, but for molecular clouds in general. 
Thus, it seems unavoidable to think that the molecular cloud 
supersonic linewidths, rather than being hydrodynamical 
turbulent motions, are what we refer to as hierarchical grav- 
itational collapse, i.e, the local gravitational contractions oc- 
curring throughout the whole cloud, which, furthermore, is 
itself collapsing. In other words, the kinetic energy gained 
during the hierarchical collapse must come from the gravi- 
tational energy released, and therefore develops a virial-like 
relation, except that the velocity dispersion is given by eq. 
((6| rather than by the Virial relation 5v 2 /R = GE. 

The assumption of gravity driving the chaotic mo- 
tions in multiple local centers of collapse, and thus de- 
veloping a pseudo- Virial state (eq. [B}, implies that mas- 
sive, compact cores should develop larger velocity disper- 
sions for larger column densities (N ~ 10 23 -10 25 cm " 2 , 
IShirlev et al1l2003l : iGibson et al1l2009l : IWu et alJlioioh . In 
Fig. [3] we show, in the velocity dispersion-size space, lines 
of constan t column den sity according to eq. ((6]), the lo- 
cus of the iLarsonl ljl98lh relation, eq. ((4|, and the region 
where the massive cores are located. From this figure we 
note that typical local clouds with mean column densities 
of 10 21 -10 22 cm will necessarily be close to Larson's rela- 



tion, as observed in the data (e.g., lLarsoniri98ll ; lHever et all 
2009). Column densities far below this relation will be un- 
able to self-shield against background UV rad iation (e.g., 
Irlartmann. Ballesteros-Paredes fe Bergir] (|200ll )). and will 
be rapidly dissociated, or do not form at all. However, as 
discussed previously, column densities far above Larson's re- 
lation do exist - in the recently massive compact cores, which 
occupy the locus 0.1 ^ r/pc ^ 1, 1 ^ 5f/km s _1 ^ 10 (see 
Figs. [T] and [3]) and, although they do not fall on Larson's 
rela tion, they do fa ll on Heyer's. 

ILarsonl (|l98ll ) argued that supersonic hydrodynam- 
ics must be important in cloud structure, and that the 
clouds cannot have forme d "by simple gravitational col - 
laps e" . The model clouds b y Vazquez-Semadeni et al.l (|2007l ) 
and iHeitsch fe Hartmannl l|2008h do indeed exhibit impor- 
tant effects of hydrodynamic turbulence, especially in early 
stages of formation. However, hierarchical gravitational col- 
lapse comes to dominate at late stages; it is just not "simple" 



4 Note that the resulting motions of such hierarchical gravita- 
tional collapse is, in a way, supersonic turbulence. However, rather 
than being the physical ingredient that opposes gravity, it is a 
consequence of gravity itself. 
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Figure 3. Velocity dispersion-size space. The lines are the loci 
of cores and clumps at constant column density if the velocity 
dispersion is of gravitational origin, such that eq. (0 is valid. 



because the geometry is complicated, with several local cen- 
ters of collapse within the parent cloud that is contracting 
as a whole, giving rise to chaotic density and velocity fields. 
The local collapse occurs because the initial turbulence, in 
combination with strong radiative losses induces nonlinear 
density fluctuatio ns that have shorter free-f all times than 
the parent cloud (|Heitsch fe Hartm ann 2008). 

The scenario of hierarchical gravitational collapse is 
similar in spirit to the classic notion of gravitational frag- 
mentation (Hoylc 1953), except that it takes place on a sub- 
strate populated with nonlinear density fluctuations pro- 
duced by the initial turbulence. This eliminates the si- 
multaneity of the collapses that would be expected in a 
uniform medium. In this case, the driving force is grav- 
ity, as opposed to previous scenarios where the confining 
force for the (hydros tatic) clouds was the pressure externa l 
to the clouds (e .g., lElmegreenl 1 19891 ; iMcKee fe Tanl 120031 ; 
iField et al.ll2010l ). It is interesting to note, however, that in 
those works the presence of such a confining pressure was as- 
sumed, and diverse mechanisms had to be devised in order 
to attain the necessary h igh pressures, su ch as the weight 
of an atomic component |Elmegreenll 19891 ) or a recoil pres- 
sure from the ph otodissociation regions around the clouds 
|Field et al.ll2010T ). In our scenario, such confining pressure 
does not exist, and the large internal pressure of the clouds 
is simply a consequence of the ongoing collapse. 



3.2 Avoiding the star formation rate conundrum 

The idea that GMCs and their substructure may be in 
a state of gravitational collapse is not new. In fact, it 
was the first propo sed explanation for the observed su- 
personic linewidths (|Goldreich fc Kwarjfl'973 '). However, as 
it is well known, it was q uickly deemed untenable by 
IZuckerman fc Palmer I (|l974l ) , who argued that if all molec- 
ular gas in the Galaxy were in a state of free-fall, the star 
formation rate (SFR) would be of order M mo i/iff, where 
Afmoi is the total molecular gas mass in the Galaxy, and tg 
is the mean free-fall time of the molecular gas. This results 
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in a star formation rate roughly 100 times larger than the 
observed value of a few Mq yr -1 . 

This simple reasoning, however, does not necessarily ap- 
ply for real molecular clouds, since it neglects their complex, 
highly fragmented nature of their density distribution. Be- 
cause of the existence of a wide distribution of local free-fall 
times in a GMC, the densest clumps collapse significantly 
sooner than the GMC at large, and begin forming stars be- 
fore the bulk of the GMC has collapsed, as observed in sim- 
ulations fe.g. | Vazquez-Semadeni et al.ll2007l; iHeitsch et alj 



120081 ; iBaneriee et alj|2009l ; IVazquez-Semadeni et al.ll2010l ) 
The feedback from the stellar products (in the form of out- 
flows, winds, ionizing radiation, and SNe), can then in- 
terrupt the star formation (SF) process before it exhausts 
the entire mass of the GMC. Furthermore, since it is well 
known that most of the molecular gas mass is in the large- 
scale structures, this implies that the SF efficiency will re- 
main small, because the stellar activity originates from the 
objects that collapse first, which contain the minority of 
the mass. Indeed, numerical simulations of GMC formation 
and evolution including self-co nsistent stellar feedback by 
IVazquez-Semadeni et all l|2010h show that the SF efficiency 
can be maintained at realistic values throughout the evolu- 
tion of the cloud, still within the context of large-scale grav- 
itational contraction. As mentioned above, the youth and 
strong coevality of stellar populations in the clouds, suggest 
that stellar feedback acts quickly to suppress SF at any given 
location, and acts to either relocate it, for a next generation 
of SF at a different location, or even to completely disperse 
a cloud. Although the details of this process still need to be 
refined, it is clear that the SF conundrum does not need to 
apply for GMCs in general. 



4 CONCLUSIONS 

In the present contribution we have discussed the nature of 
the velocity dispersion-size relationship. We showed that re- 
cent o bservational results do not follow the standard iLarsonl 
l|l98ll) velocity dispersion-size relation. Instead, (i) mas- 
sive cores occupy the locus defined by 0.1 ^ R/pc ^ 1; 
1 ^ Sv/km s" 1 ^ 10, and (ii) rather than a single 5v — R 
relation, the entire dataset for which in dependent mass es - 
timates are available seems to follow the lHever et~aH (|2009l 'l 
scaling, Sv oc T}^ 2 R}^ 2 . We showed that these results are 
consistent with molecular clouds in a process of hierarchical 
gravitational collapse, i.e., molecular clouds collapsing as a 
whole while, at the same time, their cores collapsing locally, 
creating a complex supersonic velocity pattern that is how- 
ever not fully random, as in the standard notion of turbu- 
lence, but rather contains a dominant globally contracting 
mode. As a consequence, this turbulence cannot oppose the 
contraction, but rather feeds from it. 

We emphasized that, although hydrodynamic turbu- 
lence in the warm ISM must play a role in producing the ini- 
tial m olecular cloud and its condensations (|Clark fc BonnelH 
120051 ). once gravity dominates motions at late stages natu- 
rally results in the pseudo-Virial relationships seen in the 
data. Thus, it is not necessary (let alone likely) that clouds 



and their cores are in pressure equilibrium with the external 
medium, nor is it necessary to resort to some unspecified 
mechanism that can accurately keep massive clouds with 
complex, non-spherical, and clumped density distributions 
in an approximate equilibrium for many crossing times. 
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